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ABSTRACT 

We use a one dimensional hydrodynamical code to study the evolution of spher- 
ically symmetric perturbations in the framework of Modified Newtonian Dynamics 
(MOND). The code evolves spherical gaseous shells in an expanding Universe by em- 
ploying a MOND-type relationship between the fluctuations in the density field and 
the gravitational force, g. We focus on the evolution of initial density perturbations 
of the form Si ~ r~ s for < s < 3. A shell is initially cold and remains so until 
it encounters the shock formed by the earlier collapse of shells nearer to the centre. 
During the early epochs g is sufficiently large and shells move according to Newtonian 
gravity. As the physical size of the perturbation increases with time, g gets smaller 
and the evolution eventually becomes MOND dominated. However, the density in the 
inner collapsed regions is large enough that they re-enter the Newtonian regime. The 
evolved gas temperature and density profiles tend to a universal form that is indepen- 
dent of the the slope, s, and of the initial amplitude. An analytic explanation of this 
intriguing result is offered. Over a wide range of scales, the temperature, density and 
entropy profiles in the simulations, depend on radius roughly like r 5 , r~ 15 and r 1,5 , 
respectively. We compare our results with ~XMM-Newton and Chandra observations 
of clusters. The temperature profiles of 16 observed clusters are either fiat or show 
a mild decrease at R }t 200kpc. MOND profiles show a significant increase that can- 
not reconciled with the data. Our simulated MOND clusters are substantially denser 
than the observed clusters. It remains to be seen whether these difficulties persist in 
three-dimensional hydrodynamical simulations with generic initial conditions. 

Key words: cosmology: theory, observation, dark matter, large-scale structure of 
the Universe — gravitation 



1 INTRODUCTION 

In the standard cosmological paradigm the Universe is pre- 
dominantly filled with dark matter (DM). The success of this 
paradigm is indisputable. It explains a multitude of observa- 
tional data over a wide range of scales, from rotation curves 
of local galaxies to the anisotropies of the cosmic microwave 
background (CMB) at high redshifts. 

Despite this success, several fundamental problems re- 
main to be resolved. Although the best evidence for DM is 
the flattening of rotation curves of many spiral galaxies, the 
scenario does not explain the shape of these curves in the in- 
ner regions of some galaxies (e.g de Blok et al. 2001, Hayashi 
et al. 2004). The DM scenario predicts specific correlations 
between the internal properties of clusters of galaxies (e.g. 
the mass vs temperature, entropy vs temperature). While 
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such correlations are observed, their parameters do not agree 
with the model predictions. There are various puzzles related 
to the evolution of the galaxy population as well. All these 
problems and several others can probably be attributed to 
physical processes unrelated to the dark matter. But they do 
allow us some leeway in exploring other alternatives to the 
DM scenario. One possibility is to superimpose additional 
forces acting only in the dark sector (e.g. Gradwohl & Frie- 
man J.A. 1992; Farrar & Peebles 2004; Sealfon et al. 2005; 
Shirata et al. 2005). These forces can work in the direction 
of boosting the clustering on small scales, leaving the larger 
scales intact. Several desirable features follow from this sce- 
nario (e.g. Nusser, Gubser & Peebles 2005). 

Another approach is to modify Newton's g N ~ 
l/(distance) 2 law for the force of gravity, relinquishing the 
DM scenario altogether. This approach was originally in- 
troduced to account for the flattening of rotation curves 
without invoking dark matter (Milgrom 1983; Bekenstein 



2 A. Nusser & E. Pointecouteau 



& Milgrom 1984). After all, the DM particle has yet to be 
discovered (Bertone et al. 2004) and we lack an experimen- 
tal verification of Newton's laws at low accelerations. The 
essence of Modified Newtonian Dynamics (MOND) is to re- 
place Newton's law of gravity at sufficiently low accelera- 
tions g N < g„ l by g ~ ^%g^ where g « 1.2 x 10" 8 cm s" 2 
is found to give the best results in fitting the rotation curves. 

Recently a few attempts have been made to confront 
MOND with observations of the large scale structure (LSS) 
in the Universe (McGaugh 1999, Sanders 2001, Nusser 2002, 
Knebe & Gibson 2004, McGaugh 2004). When MOND is ap- 
plied to a uniform background it predicts the collapse of any 
finite region in the Universe regardless of the mean density 
in that region (e.g., Felten 1984, Sanders 1998). To solve this 
problem Sanders (2001) proposed a two-field Lagrangian 
based theory of MOND in which the Friedmann-Robertson- 
Walker (FRW) background cosmology remains intact in the 
absence of fluctuations. He argued that this theory leads 
to LSS resembling Newtonian dynamics with CDM-like ini- 
tial conditions. More recently Bekenstein (2004) proposed a 
general relativistic version of MOND in which the behaviour 
of cosmological background is very close to a FRW. Nusser 
(2002) employed the "Jeans swindle" (Binney & Tremain 
1987) to write a MOND type relation between the fluctu- 
ations in the density and the gravitational force field. The 
relation can be derived from a Lagrangian and is equivalent 
to Sanders' two-field theory in the limit of small coupling 
(Sanders 1998). And it seems that a similar relation can 
be derived from Bekenstein's (2004) theory as well. Nusser 
(2002) then implemented this relation in a collisionless N- 
body code to simulate the evolution of LSS under MOND. 
This work showed that MOND, albeit with go smaller than 
the standard value inferred from rotation curves, produces 
LSS similar to that seen in simulations of viable variants 
of the Cold Dark Matter (CDM) scenario. Knebe & Gibson 
(2004) used high resolution simulations of collisioness parti- 
cles to study the LSS and "halo" profiles in MOND. Accord- 
ing to these authors, MOND leads to reasonable predictions 
for the LSS even for the standard value of go. Knebe & Gib- 
son also examined the properties of halos and concluded that 
density profiles in MOND have similar shape to the profiles 
seen in Newtonian simulations. 

A MONDian Universe is dominated by ordinary bary- 
onic matter made mainly of primordial hydrogen and he- 
lium. Baryons in cosmological systems have a relatively short 
mean free path, and should be treated as a hydrodynamical 
fluid. This fluid can collapse under its self-gravity, shock- 
heat during the collapse, cool, and form stars. Reliable three 
dimensional simulations of these processes are difficult and 
very CPU demanding, even in Newtonian dynamics where 
the Poisson equation can be solved efficiently. The problem 
is greatly simplified if one restricts the analysis to symmet- 
ric perturbations. Stachniewicz & Kutschera (2004) adopted 
this approach to the collapse of low mass objects at high red- 
shifts using a spherically symmetric hydrodynamical code 
with MONDian gravity. In this paper we address the forma- 
tion of clusters of galaxies in MOND using hydrodynami- 
cal simulations of spherically symmetric perturbations hav- 
ing initial power law profiles. The perturbations are evolved 
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using a one dimensional Lagrangian hydrodynamical code. 
Cooling and heat conduction are not included in this code. 
These processes are important in galaxies and small groups 
of galaxies. However, they are likely to play a minor role in 
governing the evolution of the gaseous intercluster medium 
(ICM) in the outer regions of massive clusters. Several works 
have confronted MOND with observations of clusters (e.g. 
The & White 1988; Aguirre et al. 2001). Aguirre et al. (2001) 
assumed the hydrostatic equilibrium to derive MOND tem- 
perature profiles from X-ray data of three nearby clusters. 
They concluded that MOND predictions for the tempera- 
ture profiles disagree with observations. Our approach here 
is to derive cluster properties resulting from the dynamical 
evolution of initial perturbations. This approach will enable 
us to make detailed predictions independently of the obser- 
vations, and to test several key assumptions such as that of 
the hydrostatic equilibrium. Therefore, we expect it to pro- 
vide stringent constraints on MOND from the observational 
data. 

The remainder of the paper is organised as follows. The 
notation, the equations of motion, and the initial conditions 
are described in J21 An analytic treatment of the evolution 
when the perturbations are small is given in H2.2I The nu- 
merical model is outlined in [0 In tests using known 
self-similar solutions are presented. The results for MOND 
are described Sj^] In Sjfjour results are confronted with recent 
observations. A final summary and discussion are presented 
in ffl 



2 THE MODIFIED COSMOLOGICAL 
EQUATIONS OF MOTION 

The background FRW cosmology is described by the scale 
factor a(t) normalised to unity at the present, the Hubble 
function H(t) = a/a, and the total mean background mat- 
ter density ptot — Pdm + pb , where pd m and pb are the mean 
densities of the dark and baryonic matter, respectively. In 
MOND the Universe is made of baryonic matter and so we 
take pdm = 0. Further, we assume a vanishing cosmological 
constant since otherwise the cosmic age would be unrealis- 
tically long if pb is fixed by nucleosynthesis. We also define 
fib = Pb/pc, where p c = 3H 2 /(8itG) is the critical density. 
These cosmological quantities are related by Einstein equa- 
tions of general relativity. Let r and x = r/a denote, respec- 
tively, physical and comoving coordinates. The fluctuations 
over the uniform background in the matter distribution are 
described by the comoving peculiar velocity v = dx/ dt of 
a patch of matter, the density contrast 5(x) = p(x)/p — 1, 
where p(x) is the local density, and the fluctuations in the 
gravitational force field, g. Further, the thermal state of the 
gas is described by the pressure, P and the internal energy, 
u, which are assumed to be related to the density by a per- 
fect equation of state P = (7 — l)pu, where 7 = 5/3. Here, 
we find it convenient to work with p = Pj p, instead of the 
physical pressure P. The Newtonian equations of motion 
governing the evolution of a spherical perturbation are: the 
continuity equation 

™ + l_ dx[x 2 {1 + 8)v]=0, (l) 
the Euler equation of motion, 
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dv 1 1 d x p 

— + 2Hv = -g- 
dt a 



the energy equation, 

du _ p d8 
~dt ~ (l + 8) 2 ~dt ~ 



1 + 8 5 



3Hu 



(2) 



(3) 



where the "pressure" , P, and internal energy, it, are related 
by the equation of state, 



P= ( 7 -l)(l + 5)t 



(4) 



The final equation of motion that is needed is a relation 
between g and the density, 8. In the Newtonian theory g = 
gr N where g N satisfies the Poisson equation, 



-^-d x (x 2 g N ) = ~4TvGap b S = -^-aSl h H 2 8 



(5) 



In MOND g = g M where g M is related to g N by g N = 
g M fi(g M /go) with /%) = y(I+?/ 2 ) 1/2 (Milgrom 1983, Beken- 
stein 2004). This relation yields 



1 




(61 



Neglecting thermal effects, the linear Newtonian theory 
implies that \g\ ^> go at sufficiently early times (see t|2.2H . 
But \g\ decreases with time and MOND eventually takes 
over the evolution of the perturbation. Matching MOND to 
observations of rotation curves of galaxies gives go — 1 x 
I0" 10 m s~ 2 . We further assume that go is constant with 
time. 



2.1 Initial conditions 

We chose a power law for the initial density contrast profile. 
We express the mean density contrast inside a distance x 
from the centre of symmetry as 



5i(* = 0) = <5 f— V 

\Xo/ 



(7) 



when normalised to redshift zero (z = 0) according to New- 
tonian theory. We choose <5o = 1.68 so that xo is roughly 
the Lagrangian size of the collapsed object at z = 0. This 
physical meaning for xo is only valid in Newtonian dynam- 
ics. This way of normalising the perturbation has no bearing 
on the final results for MOND and is adopted here only for 
convenience. 

The initial conditions are given at some very high red- 
shift, z\ as 



li(zi) 



0) 



(8) 



where D is the (Newtonian) linear growth factor from z\ 
until z = (e.g. Peebles 1980). The initial velocities are set 
according to the growing mode of linear theory as 



v 1 (x) = — HiX , 

where Hi = H(zi). The internal energy is 
ui(x) = . 



(9) 



(10) 



2.2 The limit of small density perturbations 

Here we treat the limit of small density perturbation when 
i < 1 such that the convective terms in the equations of 
motion can be ignored. We will show that in this limit the 
density contrast tends to the form S <x 1/x independently 
of the amplitude and slope of the initial contrast density 
8i oc x~ s . We make the following assumptions, (i) fib = 1 
as is the case if the analysis is restricted to sufficiently early 
times, (ii) the gas is initially cold and remains cold as long 
as 8 <C 1 , and (Hi) the transition to the MOND regime is 
abrupt so that g = g N ^g n /\g N \ for g < g , and g = g N 
otherwise. 

The initial mean density contrast, 8i = 8o(x/xo)~ s , is 
given at some very high redshift z\. We first assume that 
8 > 1. In this case \g\ oc x8i ~ x 1 '" decreases with radius, 
so that at redshift zo, there exists a distance x m {z) such that 
all shells with x ^ a: M (zo) are in the MOND regime, while 
those with x < x M (zo) axe still in the Newtonian regime. 
We compute below the density profile for x 2> x M . The 
reason for excluding shells at x ^ x M will become clear as 
we proceed. We assume that a shell at x 3> x M enters the 
MOND regime at z M (x) < z\. For zi > z > z M all shells 
interior to x are in the Newtonian regime and the amplitude 
of the peculiar gravitation force field at x is given by 



1 + * 



,8{z)GM 



(11) 



where 8{z) is the density contrast at z, and M — (4ir/3)px s , 
p being the background density at redshift zero. In the New- 
tonian regime, we write 8(z) = (1 + zi)8i/(l + z). Therefore, 

4-7T . 



-Gp8 1 {l + Zl ){l + z)x 



(12) 



The transition to the MOND regime occurs at z M at which 
| g | = go- Therefore, 



1 + ^ 



47T ,_ , _ X 

— (l + z^SiGp— 

■5 go 



(13) 



The density contrast in the MOND regime grows like 1/(1 + 
2) 2 for sufficiently late times (z <C z M ) (Nusser 2002). We 
can use this result if x 3> x M (zo) since this condition implies 
that zo <§C z M . The density at z M is (1 + z\)Si / (1 + z M ) , which 
must be multiplied by (1 + za/) 2 /(1 + zo) 2 to obtain the 
density, 8(zo), at zo- In this last step we have assumed that 
the mass of matter encompassed within a sphere of radius 
x M , which is Newtonian, does not contribute significantly to 
8 at x x M . The result is 



8(x,zo) = 



3po 



1 



(14) 



4-KGp (1 + zo) 2 x ' 

for x 3> x M and s > 1. For s < 1 one can show that the last 
expression applies at x <S x M . 

3 NUMERICAL MODEL 

The evolution of a spherical perturbation in an otherwise 
uniform Universe is followed using a one-dimensional La- 
grangian hydrodynamical code. The code follows the evolu- 
tion of of iV shells inside a region of initial comoving radius 
i m „. At the initial time, the shell i + 1/2 is bounded by the 
Xi — (i — 1) da; and Xi+i = Xi + dx where da; = a; max /iV. 
The shell is then assigned a "mass", 
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dM t+1/2 = y [xi+i(l + 5i(x i+1 )) - 4(1 + - (!5) 

where 5i is given by ©. This mass is maintained constant 
throughout the evolution ol the shell. The initial velocities 
at x = Xi are set according to JjJJl. A leapfrog integration 
scheme is used and so the positions, x™ , at the time step t = 
r and the velocities, v n+1/2 , at t = t n ~ 1/2 = i™— dt n ~ 1/2 /2, 
are advanced as follows, 



-1/2 



-1/2 



1 + H n dt H 

MX?) 2 P?+l/2 -P7-1/2 _ £_ 



(a") 2 



dt" 



n + l n i n+1/2 j ,n 

= + V, at 



where 



at" = i(dr +1/2 



-1/2 



-1/2n 



(16) 
(17) 

(18) 
(19) 



dMi = - (dMi- 1/2 + iMj +1/2 ) , 

and p = p + q, where q accounts for artificial viscosity. Cast- 
ing the standard expression for q (e.g. Richtmyer & Morton 
1967, Thoul & Weinberg 1995) in terms of comoving coor- 
dinates gives, 

9"+ + i/2 = 2cq(H^ ) 

1+1/2 (i+<5r + i 1 / 2 )- i + (i+^ 1/2 )- 1 

\ V ^ 2 -vr i/2 \{v^' 2 -vr i/2 ), (2i) 

v n+i/2 < ^ and q n+i^ = g otnerwise . H ere 



if ^»+l/2 



n + 1/2 „?i + l/2 

- a — (^ri 1 + x?) + ^-^ v ? +i/2 



fjn + l/2"i (^) 

is equal to the actual total physical velocity and we take 
c q — 4 which would spread the shock over ~ 4 shells. In the 
Newtonian case the gravitational force field is 



9i = 



8 Ti 



(23) 



where the summation is over all j having Xj < Xi. This 
expression is used in @ to obtain the gravitational force 
field in MOND. 

The density and energy are updated, respectively, ac- 
cording to 

3dM i+1/2 



1 -l- 

1 + °*+l/2 



and 



n + l 
•^+1/2 



^ [(^+ + i ) 3 - (^r +1 ) 3 ] ' 



„ 2-m n+l/2 dt n+1/2 

Wi+l/2 2 + 3J/ „ + l/2 dt«+l/2 
(p" +1 +p n ) 



(24) 



+ <7i+l/2 
1 



(25) 



Cd , 



V 1+<s i+i/a 1+<5 r + i/ 2 ; ' 

The time-step dt n+1//2 is chosen to be the minimum of 

(26) 



cc 



Xi-l 



\/7(7- !)«» ' 



and 



Si-i 



«i-i 



(27) 



(28) 



where we take c d = 0.001, c c = 0.2 and c„ = 0.05 (Thoul & 
Weinberg 1995). 

The boundary conditions are: Vi = i — 0, and ap K+1/2 = 
0. We have not included any explicit force softening in the 
code. For the number of shells we considered here the runs 
in all cases were completed within a few days of CPU time 
on a G4 machine. 



4 TESTS OF THE NUMERICAL SCHEME IN 
THE NEWTONIAN CASE 

Since the change from Newtonian gravity to MOND is trivial 
(see Eq.0, it suffices to confront the code with self-similar 
solutions in the Newtonian case. For an adiabatic gas, the 
initial conditions given in Eq. lead to self-similar evo- 

lution only in a flat Universe with Qb — 1. Therefore, in this 
section we adopt this value of fib- Here show results from 
the simulation output at redshift, z = 0.1. Instead of scales 
variables (e.g. temperature and distance relative to their re- 
spective virial values) we work with physical quantities. The 
temperatures are given in degrees Kelvin, number densities 
in cm -3 and distances in kpc. The advantage of using physi- 
cal quantities is that they provide a more direct comparison 
with MOND where the mean density within the virial radius 
is not necessarily proportional to the background density. 

We have run the code with xo — 10 Mpc for three values 
of s: s = 0.25, 2 and 2.8. We use N = 500 shells for the runs 
with s = 0.25 and s = 2, and N = 100 shells for s = 2.8. The 
reason for the smaller number of shells for s = 2.8 is that 
the code is substantially slower in this case as the time steps 
become extremely short at the later stages of the evolution. 
All simulations were run with an adiabatic index of 7 = 5/3. 

In Newtonian dynamics, given the initial conditions 
Eq. I8I10H . the evolution of a shell can be described qualita- 
tively as follows. In the early stages the shell expands until 
it reaches its turnaround radius at which the total physical 
velocity, aHx + av, is zero. The shell then collapses towards 
the centre until it encounters the shock, transforming most 
of its kinetic energy into heat and become part of the hot 
region. 

In Fig.^we plot the temperature (top panel) and den- 
sity (bottom) versus the physical distance, 7?, from the out- 
put of the simulations at redshift z — 0.1. The shock po- 
sition, Rsh, in the various cases corresponds to the abrupt 
change in the density and temperature. The transition in 
the shock is smoother for s = 2.8 (dashed curve) as a re- 
sult of the smaller number of shells in this run. The density 
variation across the shock for s — 0.25 and s — 2 is close to 
(7 + l)/(7 — 1) = 4, as expected in strong shocks. Because 
of the smaller number in the s = 2.8 case, the density jump 
at the shock is weaker in this case. 

At each time, t, there is a shell that is currently at its 
turnaround radius. This radius is (in physical units), 

i? ta = zo(l + z)- (1+1/s) (j5o) [£) ■ (29) 
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Figure 1. The temperature (top) and density (bottom) profiles 
at redshift 2 = 0, obtained from the code for a power law initial 
profiles (see Eq[SJ for xq = 10 Mpc. The various curves corre- 
spond to different values of the power law index, s, as indicated 
in the figure. The runs are for flat Universe with Q b = 1 and no 
dark matter. All curves agree well with the self-similar solutions 
of Chuzhoy & Nusser (2000). 



We have inspected the output of the simulations at various 
times and confirm that the turaround radii agree very well 
with For s = 0.25 and s = 2, Chuzhoy & Nusser (2000) 
(hereafter CN2000) provided full self-similar solutions (see 
their figure 2). The shapes of our profiles match extremely 
well the solutions of CN2000. We also checked if the shock 
radius relative to the turnaround radius agrees with the self- 
similar solutions. At redshift z — 0.1 we find, from 12911 . 
that R u = 6.155 Mpc and 21.9 Mpc for s = 2 and 0.25, 
respectively. Inferring the shock radii, Rsh, from Fig. Q we 
get ir sh /i? ta « 1.8/6.155 = 0.29 for s = 2 and 0.57/21.9 = 
0.026 for s = 0.25. These ratios agree very well with those 
inferred from figure 2 of CN2000. 

For s = 2.8 and s = 2, the asymptotic logarithmic 
slopes near the origin agree very well with the analytic slopes 
presented in Table 1 of CN2000. For s = 0.25, neither our 
curves nor those in figure 2 of CN2000 behave according 
to the asymptotic exponents in that Table. This is because 
the analytic exponents in the s = 0.25 case are realized at 
extremely small fractions of the shock radius (see CN2000) . 



cosmological parameters of choice, the age of the Universe 
« 14Gyr. 

The initial redshifts in all the runs is chosen sufficiently 
high so that all shells in the simulations are still Newtonian. 
We investigated various values for the parameters s and xo 
(see Eq.|SJ. In Fig. Inland Fig.[3]we show the temperature and 
the local baryon number density (in cm -3 ) versus the phys- 
ical distance from the centre, R, for three different values 
of s, but the same xo- The normalisation of the initial den- 
sity profile is such that the average density contrast within 
a comoving distance xo is the same for all three values of s. 
There is a striking similarity in the shapes of the tempera- 
ture and density profiles, respectively. This is in agreement 
with the analytic considerations in section § 12.21 The pro- 
files are also insensitive to xo as seen in Fig. [3] which shows 
results for s = 1 but different xo- Even a change of a factor 
of 30 in xo fails to amount to any significant variation in the 
shape. The overall amplitude of both the temperature and 
the density profiles is also insensitive to aio. The density and 
the temperature show, respectively, ~ R~ 15 and ~ R 5 , 
dependence on radius. 

The "entropy", S = kT/n 2/3 , is plotted in Fig. H The 
entropy profiles corresponding to all the simulation outputs 
presented in the previous two figures are shown. In the shock 
heated regions the curves are practically indistinguishable at 
radii R £ 200kpc. At these radii, S cc r 1-S . 

The redshift evolution for the case s — 1 and xo = 
1 Mpc is explored in Fig. [S] showing profiles at three differ- 
ent redshifts. There is an increase by a factor of about 1.8 
in the shock location between z = 0.5 and z = 0.1, but the 
profiles at these two redshifts have similar shapes. Over the 
distance scales in the plot the temperature shows a signif- 
icant change from z = 3 until the present. But the overall 
shape at z = 3 is similar to the outer portion of the low 
redshift profiles. This indicates that the shape is preserved 
with redshift. The density profiles inside the shocked region 
at the three redshifts are remarkably similar, indicating that 
the system is in hydrostatic equilibrium. 

At early epochs near the initial times, perturbations 
evolves according to Newtonian dynamics and then, as \g\ 
decreases, they succumb to MOND. This is as long as the 
perturbations are small. Fig. |S] demonstrates that in col- 
lapsed regions the perturbations become Newtonian again. 
The figure shows \g\ relative to go as a function of R. In all 
cases \g/go\ > 1 in the inner regions. 



5 RESULTS WITH MOND 

All results here will be presented for go — lx 10 - °m s -2 and 
Qb = 0.044. This value of is consistent with the obser- 
vations of the deuterium abundance relative to hydrogen in 
Lya clouds at high redshift (Kirkman et al. 2000). It is also 
consistent with the ratio of the first to second peak heights 
in the angular power spectrum of the CMB anisotropies 
(e.g. Spergel et al. 2003). We further take h = 0.7 for the 
present Hubble constant (in units of 100 km s _1 Mpc _1 ). We 
work with a vanishing cosmological constant since otherwise 
the age of the Universe will be unreasonable large. For the 



The insensitivity of the profiles to the initial conditions 
is consistent with our findings in H2.2I A perturbation en- 
ters the MOND regime when the density fluctuations are 
still small. During this stage of the evolution they acquire 
the form 1141 independently of the initial density profile, as 
shown in H2.2I Later evolution, when the density fluctua- 
tions are large, modified the form given in H2.2I but can not 
restore the dependence on the initial conditions. It is inter- 
esting that the evolved form of the density profile in MOND 
is ~ r which is the same as the one obtained in the self- 
similar collapse of an initial 1/x perturbation in Newtonian 
dynamics (e.g. CN2000). 
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Figure 2. The temperature (top) and the density (bottom) pro- 
files in MOND for s = 0.5, 1 and 2.8, as indicated in the figure. 
The normalisation is such that the mean value of the density 
within a comoving distance xq = IMpc is the same at the initial 
time for all three values of s. 



MOND 




R (kpo) 

Figure 3. Profiles in MOND for s = 1, but for three values of 
the normalisation as indicated in the figure by the parameter xq. 



6 THE OBSERVATIONAL SITUATION 

We confront the temperature and density profiles of MOND 
with recent X-ray observations of galaxies clusters. We 
use the recently published temperature profiles by Pointe- 
couteau et al. (2005) obtained for a sample of 10 nearby 
(z < 0.15) clusters with the XMM-Newton satellite. To give 
a fair idea of the variation in shape of the observed tem- 
perature profiles, we also add the temperature profiles of 
eleven nearby (z < 0.23) clusters obtained by Vikhlinin et al. 
(2004) with the Chandra satellite. Both samples cover a wide 
range of temperatures from 2 to 9 keV. We scale each tem- 
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Figure 4. Entropy profiles in MOND at 2 = 0.1 for various 
values of s and xq. The slope, d In S/d In il, at R ^ 10 3 kpc is 
1.3, approximately. 
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Figure 5. Profiles in MOND for s = 1 and xq from outputs at 
three different redshifts, z, as indicated in the figure. 



perature profile to its value at 500 kpc (i.e kT(r — 500 kpc)) 
in order to compare them with our simulated profiles. To 
avoid any extrapolation, we kept only the 16 clusters which 
were observed beyond 500 kpc (i.e 8 XMM-Newton and 
8 Chandra clusters). The scaled temperature profiles are 
shown in Fig. |7| (top panel) . 

We also compare the observed density profiles with 
MOND's profiles. The observed profiles correspond to the 
parameterised form derived by fitting the X-ray surface 
brightness profile of each cluster. For this comparison, we 
only use the XMM-Newton profiles. The best fit models for 
the density profiles of the ten clusters from Pointecouteau 
et al. (2005) are plotted together with the MOND density 
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Figure 6. The ratio of the absolute value of the peculiar gravita- 
tional field relative to go = 1 X 10 — 10 ms~ 2 as a function of radius 
for runs with the same normalisation and three different values 
of s. 



profile in Fig. [7| (bottom panel). All profiles are scaled to 
their respective values at 500 kpc. 

All of the observed clusters are identified as relaxed 
and show small deviations from sphericity. The temperature 
drop towards the centre in the observed profiles is the result 
of cooling. Since cooling is not included in our simulations, 
we restrict the comparison to R ^ 200 kpc. 

For better visualisation, all profiles are scaled to their 
respective values at R = 500 kpc. At that radius, the mini- 
mum, maximum, average, and median values of the observed 
temperature are (in 10 7 K), 1.9, 11.9, 6.6, and 6.8, respec- 
tively. The corresponding values for the baryon number den- 
sity are (in 10~ 4 cm~ 3 ), 4.4, 20.6, 13.2, and 18.2, respectively. 
At this radius MOND yields a temperature range consistent 
with the observations. However, among the MOND runs of 
all values of s and xo that we shown, the minimum value 
number density at 500 kpc is MOND density in MOND is 
0.07 cm -3 . This is substantially higher that the maximum 
observed number density of 0.0026 cm~ 3 . More importantly, 
the temperature profile, ~ r ' 5 , in MOND is, in clear con- 
trast to the observed profiles, either flatten or show a mild 
decline in the outer regions. The entropy profile in MOND 
has a logarithmic slope of ~ 1.5, as compared to 0.94 ± 0.14 
and 0.95 ± 0.02 found in the observations (Pratt & Arnaud 
2005 and Piffaretti ct al. 2005, respectively). This is another 
manifestation of the mismatch between temperature profiles 
in MOND and the observations. 
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Figure 7. Top: MOND simulated temperature profiles (thick 
black line) compared to the individual profiles of clusters observed 
with XMM-Newton (red dashed lines) and Chandra (blue dot- 
ted lines). Bottom panel: MOND simulated density profiles (thick 
black line) shown with respect to the observed density profiles of 
nearby clusters by the XMM- Newton satellite (red dashed lines). 
In both panels the profiles have been normalised according to 
their respective values at r = 500 kpc as marked radius by dotted 
vertical line. 



7 SUMMARY AND DISCUSSION 

In this paper we aimed at: (1) studying the density and 
temperature profiles of spherically symmetric perturbations 
evolved under MONDian gravity, and (2) comparing the pre- 
dicted profiles with recent X-ray observations of the intra- 
cluster medium (ICM) . A key finding here is that the profiles 
in MOND acquire a nearly universal shape independently of 
the initial fluctuations. This is in contrast to the Newtonian 
theory where the profiles are sensitive to the shape of the 
initial profiles. The shock radius in MOND is nearly inde- 
pendent of the amplitude of the initial perturbation while in 
the Newtonian case it is proportional to the amplitude (see 
Eq.HU. 

The insensitivity to the initial density profile, leaves 
very little leeway in matching MOND predictions to ob- 
served properties of the ICM. In Newtonian dynamics, one 
can, to a great extent, resort to varying initial conditions to 
try to match the observations. In MOND one has to depend, 
almost entirely, on non-gravitational physical processes. One 
such process may be thermal conduction. The thermal con- 
duction time scale is estimated as R c ~ (2nTt /5nfcT) 1//2 . 
Assuming a thermal conductivity, k, of 0.2 times the Spitzer 
value k Sp = 9.210" 7 T 5/2 erg s" 1 K _1 cm" 2 (Spitzer 1962), 
we obtained R c ~ 200 kpc for t = 14 Gyr, n = 0.001cm" 3 
(see Fig. Hat R = 10 3 kpc), and T = 5 x 10 7 K. There- 
fore, thermal conduction can help, especially for the high- 
est temperature clusters of T ~ 10 8 K. For T < 5 x 10 7 K 
conduction seems inefficient at producing the desired tem- 
perature profile. Note that we have used the cosmic time 
t = tu ~ 14 Gyr to estimate R c . In the outer regions one 
should actually take £h — tcoiiapsc where t co ii apso is the time 
at which these regions collapsed and shock-heated. Accord- 
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ing to Fig. |^| the shock radius have grown from 600 kpc at 
z = 3 to 7000 kpc at z = 0. The time interval between z = 3 
and z = is £h — fcoiiapse = 0.25fa which reduces R c by a 
factor of 2. Non-gravitational heating of the ICM can also 
be invoked to alter the profiles, but it is hard to imagine a 
realistic mechanism that will flatten the temperature in the 
outer regions (e.g. Borgani et al. 2005). 

Our findings as well as those of previous authors, pause 
some puzzles for MOND to overcome. While Newtonian dy- 
namics has its own problems in matching the cluster data, 
the lack of dependence on the initial conditions, makes 
MOND harder to reconcile with these data. 
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